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Abstract 

Despite measurements which date more than 20 years ago, no straight- 
forward solution of the ratio of the parity-conserving (P-wave) to parity- 
violating (S-wave) decays of the hyperons has been obtained. Here we use 
two 2-point methods in QCD sum rules to examine the problem. We find 
that resonance contributions are needed to fit the data, similar to a chiral 
perturbation theory treatment. 



PACS Indices: 12.15.Ji, 13.30.Eg, 12.38.Lg, 11.50.Li 



1 Introduction 



The nonleptonic decays of the hyperons occur with pion emission, e.g., A — > 
p + n~ . Measurements of the decay rates and the S/P (pv/pc) ratios of the 
emitted pions were carried out over 20 years ago [Q. They remain of interest 
today because no one has been able to provide a relatively simple explanation 
of the S/P ratios. 

To-date a variety of approaches have been used. Some of the early work 
used a soft pion approach In this limit the S + — > nii + decay with an 
S-wave pion vanishes. The £~ — > nn^ decay amplitude can be obtained 
approximately by an adjustment of the SU(3) F/D ratio. The soft pion 
approach in the S-wave and poles in the P- waves approach (see Fig.l) was 
used by Donoghue et al.[JIJ , who argue that there could also be a direct 
coupling, as shown in Fig. Id, but they too have difficulty in fitting the S/P 
ratios. Other work is that of ref. ||. Most recently, Barasoy and Holstein 
have used chiral perturbation theory, but have had to include (70, 1~)§ 
resonances and parameters to obtain a reasonable fit to the data. 

Of the seven decays, S + — > pir , S + — ► mr + , £~ — ► mr~ , A — ► 
mr°, A — ► pn~ , H~ — > A°7r~, S° — > A°7r°, there are only four indepen- 
dent ones if isospin symmetry holds. Experimentally, the SU(3) 27-plet is 
smaller than the octet by a factor of approximately 20. Like those before us, 
we choose the 4 independent decays as those with a charged pion, namely 



K ■ 


£+ - 


-> n-K 


EI : 


£~ - 




A : 


A - 


^P7l~, 









In the present work we use the method of QCD sum rules with two 2- 
point formulations for the three-point correlators needed to obtain coupling 
constants, which we discuss in the next section. We find that in order to find 
stable solutions for the sum rules we must explicitly introduce single-pole 
resonance contributons, analogous to the addition of resonance contributions 
in Ref. @. Since this introduces new constants which can, however, be used 
to fit the data, we did not proceed to investigate the last two decays (A°_ 
and Si). 
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Fig. 1 Sigma decay into a nucleon and a pion via soft pion approach. 



2 Methodology 
2.1 QCD Sum Rules 

QCD sum rules were introduced by Shifman, Vainshtein, and Zhakarov M. 
It is a useful method to obtain properties of hadrons and it uses QCD explic- 
itly. The short range perturbative QCD is extended by an operator product 
expansion (OPE) of the correlator, giving a series in inverse powers of the 
squared momentum with Wilson coefficients. The convergence at low mo- 
mentum is improved by using a type of Laplace transform, called the Borel 
transform. The coefficients involve universal quark and gluon condensates. 
This quark-based calculation of a given correlator is equated to the same 
correlator obtained via a dispersion relation, giving sum rules from which a 
property can be estimated. The method can be extended for quantities in 
an external field, such as the magnetic coupling to a nucleon in an electro- 
magnetic field ||. 

The method begins with a correlator 

U(p) = ijd i x e iq - x < | T[r](x)f](0)} | > , (2) 

where r] has the quantum numbers of the hadron being studied. For a proton, 
we may take 



tj(x) = e a6c [n aT (x)C 7 y(a;)]7 5 7^ c (x) , (3) 
fj = e abc [u b lu Cu aT }d c Yl 5 , 
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where a, b, c are color indices and the notation of Bjorken and Drell is used. 
The quark field operators d, u, s destroy these quarks, C stands for charge 
conjugation and T for transpose. 

The "currents" rj are not unique||, but the form given in Eq. (|3]) has 
been used by many authors. The correlator can be written as an operator 
product expansion 



U = CJ + Y / C n (p 2 )O n (4) 

n 

where the operators O n can be ordered by dimension and the corresponding 
Wilson coefficients decrease by increasing powers of p 2 . 

The correlators II have structure functions IP, each of which satisfies a 
dispersion relation (P 2 = —p 2 ) 



W(P 2 ) = 1 f ImU ^ 2 ds (5) 



7T JO 



Subtraction terms in Eq.(4) are eliminated by means of a Borel transform, 
which guarantees convergence, 



B[F{p 2 )} = Km (V) (n+1) (^7— T] TF(P 2 ) (6) 

-p 2 ^oo , -p2/n^M% { P > 

1 /-oo „ 

= - ds ImF{s)e- s/M B . 

7T JO 

There are a number of ways to use these sum rules: 

(i) a two-point method with or without an external field; 

(ii) a three-point method with couplings and momentum transfers con- 
sidered explicitly. This method has fewer susceptibilities but it is more com- 
plicated; it may require non-local condensates. 

In this article we will use only two-point methods. We will compare the 
two point method in an external field with the two point method with a pion 
creation matrix element. 
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2.2 E_ in an External Pion Field 

The calculation of the non-leptonic decays of hyperons is similar to that of 
the weak pion-nucleon coupling constant. If we neglect the mass difference 
in the baryon octet, then E — > Nn is quite akin to N w ^ k Nit. The primary 
difference is that the latter is due to weak neutral currents and the former 
due to charged currents. 
We use the operators 



V x- = e abc [d aT C lil d b ]^^s c , (7) 
Vn = e abc [d b lt/ Cd aT }u c Yl 5 

for the S~ and for the n. In addition, we need the weak interaction, for 
which we use the local one, 

= _ 75 ) s s i n Q c _|_ uY(l - >y 5 )d cos6 c , (8) 

where Gf is the Fermi coupling constant and 9c is the Cabibbo angle. The 
correlator is 

Tl = iJ d 4 x e ip - x < | T[rft(x) H w f} n (x)\ \ > 7r , (9) 

where the quarks propagate in an external field. Since the point at which the 
external pion field is at zero momentum transfer 0, the three-point function 
for the vertex is reduced to a two-point function. To the order considered 
here, the QCD diagrams which contribute to the correlator are shown in 
Fig. 2. In the diagrams, the wavy line represents a W boson, the dashed 
line represents a pion. The diagrams are evaluated in momentum space. Fig. 
(2a) gives no contribution. Diagram (2e) is quite different from the others. 
It involves the weak matrix element 

< tT I J a | 7T° >= V2F^q a , (10) 

where F n is the weak pion form factor and q is the momentum of the 7r°. 
The contribution of this diagram cannot be neglected. There are additional 
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higher order terms, e.g., gluonic corrections, which we omit here. Reasonjs 
will become apparent further on. 

We obtain the following results, using dimensional regularization in 4 — e 
dimensions 

^ , 16 A < qq > p 2 Inn 2 r . . 43 e , 9 . . 15 e , «... 

Fig.2b : ™ * y [/(l - 7 e + -e -- In p 2 ) + 4m(l - ye + -e -- In P 2 )\ 

(l-7s), 

4A < qq > p 2 \np 2 r ,, 8 e , 9n , 11 e , , S1 

F^.2c : ^-[i>(l- 7 e + -e--lnp 2 )+4m(l-7e+ T e--lnp 2 )] 

(1 " 7s) 

4A < qq > jo 2 ln» 2 5 , 13 , 355 13e , 9X1/ 

Fio.2d : ^— f — - 6 ml- 7 e + e lnp 2 (1 - 75) , 11 

y (4tt) 4 l 9 ^ 4e V f 156 8 1 ;JV ' ; ' V ' 

j?- o v^i^P^lnp 2 2 6m ^ 7 e 9 2 

F ^- 2e : (4^6 [P +^( 1 + 2^ e+ 6"8 e +7s) ' 

where 7 is the Euler constant, m the strange quark mass, and 

A = V2G F sin 9 C cos fl c . (12) 

From these equations, it is clear that, as in determining the weak pion nucleon 
coupling constant ]7|, we need to include vertex renormalizations. There are 
several of these, shown in Fig. 3. For Fig. 3a we obtain 

77 = ur v Csy»y 5 , (13) 



(a) (b) (c) 

Fig. 3 Vertex renormalization diagrams 



For Fig 3b we get 



IU3&) = ~{1 ~ e -)< qq > | 7 ,(1 " 75) • (14) 
For Fig. 3c we obtain 

v = d T C^Y v u, (15) 
^ / \ 8Am , 7e e , ., 9 , c/9 . 

r " (c) = -(S)^ (1 -Y + 2 )( * 2 > - e/2 (l-75). 

Employing these 3 vertex renormalizations and calculating the diagrams of 
Fig. 2 with them, we obtain 

n = — ^-{[- * + T rn] < qq > (1 - 7fi ) - [v^^p + ..•](!+ 7 5 )} -(16) 

We have not carried out the renormalization for Fig. 2c because we shall see 
that it is not needed. 

For the phenomenological (or so-called right-hand) side we have 

11 = -AjvA s - — (A s + A P -f 5 ) 



+ jhr N ^ + ^ ) j^w^ (17) 

plus the continuum. The inclusion of resonances is indicated by the (ci+c 2 7 5 ) 
term; the mass M* represents the resonance mass. The resonance terms can 
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be separated into two parts, each of which is a single pole term; we use 
for these single pole terms. The choice of M is arbitrary; any other 



M(j/-M) 

choice will simply lead to different values for c~ and . The single-pole 
quantities c~,c^~ are functions of momentum, or after the Borel transform 
they are called c~,c^ and are functions of the Borel mass, Mb- We define 
M = \{M N + M E ) and AM = M E - M N . To first order in AM we find 

n = - {p 2 X H%y {Mp 2 + M 2 ) - A Pl5 (p 2 - M 2 ) + A S 2M £> + A P A M i> 7 e 

, x x (i> + M)(c-+c b " 75 ) 
+AwAs M(p2 _ M2) 

We abandon the /p and ^75 sum rules because Ap A M would vanish in 
the analogous pion-nucleon vertex sum rule; this term is too sensitive to 
AM. The single pole term within the square bracket is similar to that which 
occurs in determining the pion-nucleon coupling constant. Fig. 2e requires no 
renormalization for the p 2 term; we have not carried out the renormalization 
for the m fa and m ^675 terms. After inclusion of the continuum and carrying 
out a Borel transform we obtain 

A 7fi A/f 8 F 

n = -i^^[y^^ 2 (l-75) + v / 2^F.(l + 75)] 

2M 2 -K- 
= -XnX^M^T - 1) + A p^ + ( c a + c fe "75)]e m b , (19) 

where E n represents the continuum contribution, 

K = l-(l + x+- + ... + — )e"% (20) 
2 n! 

with x = s/M|, where s is the continuum threshold. For the nucleon it was 
found that s w 2.3Ge1/ 2 and for the £' s w 3.2G e y 2 §; here we usually 
take an intermediate value of s = 2.8GeV 2 for the transition, but explore 
other values for stability. Ajv and As are known from previous studies @: 
AjvAs = (27r) 4 A n A s = 0.303GeV 6 

A crucial point in the analysis of the sum rules to to recognize that 
there is some uncertainty in the p-dependence of the single-pole terms. The 



assumption used in Eq.(19) is the simplest possible. There almost certainly 
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is other M B dependence of the pole term, which makes no difference at the 
point Mb = M, but which can give stable solutions. We find the reasonable 
but not obvious behavior that in every case the single-pole constants are 
linear functons of Mb- 

As for the pion-nucleon weak coupling constant, we gain some stability 
by multiplying both sides of the QCD sum rules by M B and carrying out 
(1 — M B d/dM B )M B IT. This removes the great sensitivity in A s to M B from 
the double pole on the right-hand side, although some remains. We obtain 



-76 



M%amE 2 + y^E A M B °F n ) - ~X N ~X s c~ M 2 e~ M2 / M B r 



A _ 4 3 ?r 2 L 4 / 9 V 3 —a f 1 - n - » / --iy^-ffl--- - rtl") 

3 ~" A w A s [2(l-f^) + l]M2 e -* 2 /Ml - " 

^^{^M%amE 2 - V2E 4 M B °F n ) + \ N \^ M 2 e~ M2 / M \ ^ 

The parameters that used in this work are A = 3.57 10~ 6 GeV -2 , m=.15 
GeV, a=.55 GeV 3 , and L = 0.621 ln(10 Mg). The experimental values are 
A s = (4.27±0.02) x 10~ 7 , A P = (-1.52 ±0.16) x 10~ 7 , A S /A P ss -2.8. The 
single-pole parameters which give stable solutions for A s are 
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c - = (-6.37+10.7 M B )xlO- 7 , and for A p are = (1.37 - .25 M B )xlO- 7 . 
The solutions are shown in Fig. 4. 



2.3 S| Decay in an External Field 

For £ + we take 

v + = e abc [w aT C7XhVs c • (23) 
The relevant diagrams are shown In Fig. 5 
Fig.ha : 

2A <qq>p 2 hip 2 4 25 
^ : (4^ (3^+^(1-75) 

2A<gg>p 2 lnp 2 4 3 
^•5c : (- p+ -m)(l -75) 

F«0.5d : 

F^.5e : (24) 

„. rx 8Ap 4 lnp 2 2/i 3 25 9 . 2 , , y/1 3 73 9 , 2X1/1 , 
Fzo.5f : / xc \p (1 ~ -7e H e chip ) — 4m ^(1 76 H e - -e In p 2 ) (1 + 75) 

77 • , 8A<gg>p 2 lnp 2 3 e 2 e e 2w , 

Fi0.5# : — — [2 p(l - >ye + -e - - \np z ) + m(l - + - - -lnp z )J(l - 75) 



(47r) 4 e ' 2 2 ^ ' y 1 2 2 

qq>p 2 \np 2 5 e , 2 \ 3 e 

-7^^ [^(1 " 7e + ' 2 ln P ) + 2m(l - 7 e + -e - - 



r,. n < gg > p 2 lnp 2 5 e n 2 . rt 3 e, .-, , , 

Fm.5/i : ^-[^(l-7e+ 7 e--lnp 2 ) + 2m(l-7 e + - e --lnp 2 )](l-7 5 ) 

Fig.5i : 

. • 4A < ss > p 2 lnp 2 r 9 e. ,,,,, x 

Fig-5j : ^ ~ -76 + ^6-2 111 P )](! + 7s) 

Once again we need to carry out vertex renormalizations for Figs. 5f, g, h 
and j. We omit the details and simply show the results. As for the S~ decays 
we omit the fa and farys terms. We find 

10 Am < qq > p 2 \np 2 8 7 Ap 6 \np 2 

n = -y ^ ( 1 -^ + 8 (3-2 ) l4^ (1 +75) - (25) 

After taking a Borel transform, adding anomalous dimensions and con- 
tinuum contributions, this becomes 
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(a) 



(b) 



(c) 



(d) 




Fig. 5 Processes for E + decay into a neutron with an external 7t + field 



= W Am<qq> _ 4 4/9 _ 8 _ 7 _4 
3 (4vr) 4 v v 3 2 ; (4tt) 



:(l+ 75 )MlE 4 L- 4 / 9 .(26) 



'20 /1 t6 



A, 



A, 



43^2^4/9 v 3 



M%amE 2 + 128(1 



A JV A s [2(l-||) + l]M 2 e 
M%amE 2 - 128(1 - |^E 4 M^°) + \ N \^M 2 e 



■7) 



A /20 
43^4/9 I 3 



■M 2 /M 2 E 



Experimentally, A s = (0.13±0.02)xl0- 7 , A P = (44.4±0.16)xl0~ 7 , A S /A P 
0.003. The single-pole parameters which give stable solutions for A s are c+ 
= (-1.42 + 9.8 M B )xlO- 7 , and for A v are = (-46.3 +9.5 M B )xlO- y .The 
solutions are shown in Fig. 6. 

2.4 £ in the Pion Matrix Method 



The pion matrix method || [TTJ, [Ll| is simpler than the external field one and 
does not require any renormalizations. Instead of treating quarks propagat- 
ing in an external field with a correlator defined between the vacuum states, 
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45.0 



44. 



1.1 



1.2 



(b) 



1.3 



1.4 



Fig. 6 E + decay (a) A s , (b) A p 



the correlator is defined with a one-pion final state. This is also the start- 
ing point for light-cone sum rules which have been used for the pion form 



factor [|T2j and the pion wave function (TBI]. One makes use of the correlator 



n 



d 4 xe iq - x < vr(j9 = 0) | T[ Va {x)H w fj b {0)} | > 



(29) 



rather than Eq. (9) for the external field. For the S~ the corresponding non- 
vanishing diagrams are shown in Fig. 7. Carrying out the required algebra 
and integrations, and doing a Borel transform, we obtain 

Figla : ^A^^E 1 M 2 B < ix~ | dij s u | > ( / + 2m) (1 - 7s ) , 

Fig.lb : -^A^^E 1 Ml<n-\dn 5 u\0>(ft + 2m)(l- l5 ), 

Fig.lc : -^AF v j±y < 7i- | di l5 u | > (4mM^ 2 L" 4 / 9 - 2M 6 B E 3 L-^ 9 )(1 + 7s ) . 

Note that the contributions of Figs 7a and 7b cancel exactly, so that only 
that from Fig 7c remains and it is seen that, in contrast to the external field 



(30) 
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(a) (b) (c) 

Fig. 7 Processes for £ decay into a neutron and a 7t" by matrix method 



method, both S and P wave amplitudes have the same sign. The phenomeno- 
logic side is the same as before so that we obtain 



A a M%E Z 7 \ - A -f2 e -M 2 /M 2 „ 

A^A E [2(l-|i) + l]Me-* 2 /Mi > ^ J 



A a M%E 3 7 7 iT 4 2 p -M 2 /Ml 



where we have put F n = 1 and used the soft pion matrix element 

/-< qq > \/2a 

< o I m lsd | »- >= ^^f- = -p-^, (33) 



with f w = .093 GeV in Eqs.(0j32|). 

The single-pole parameters which give stable solutions for A s are c~ = 
( -4.3 + 12.0 M B ) x10~ 7 , and for A p are = (1.13 +3.67 M B )xlO~ 7 . The 
results are similar to those shown in Fig.4. 
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2.5 S+ Decay in the Pion Matrix Method 

The contributing diagrams are shown in Fig. 8 

FigBa : ^A^^E 1 M 2 B < n + \ ui lb d | > {£ + 2m) (1 - 75) , 

Fig.8b : - -A^^E 1 M 2 < tt+ | ui l5 d | > (2 + m)(l - 75) , 
o (47T) 

F^.8c : -4-i i <7r + |m7 5 rf|0>(^M^3 + mM^ 2 i>{(i-7 + ^ 

+ [lnMl + (l- 7 )])}(l + 75). 
F^.8d : - -A^-^E 1 M 2 < n + | ui lb d | > (2 ]£> + m)(l - 75) , (34) 
4 < ss > 

F^.8e : -A E 1 M B < tt+ | ui-y 5 d | > {£> + 2m) (1 - 75) , 

As usual, we discard the odd sum rules; in that case no renormalization is 
required. For the even sum rule we obtain 

4 Atti 

n = - — — < vr + I ui lb d I > M 2 B E x L-^{< qq > +2 < ss >)(1 - 7 5 )(35) 



3 (4tt) 



1 g A 

T < ° 1 dtl5U 1 7F+ > M ^£- 4/9 (l + 7s)(36) 
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For < ss > we take 0.8 < qq >. We thus obtain for As and Ap 




>[?fmaE l + M B E 3 ] + A^A s c+M 2 e - M2 / M i 
A]vA s [2(l-g-) + l]M2 e -MVMl 



(37) 




- c. 



.+ 

'6 ' 



(38) 



A]vA s M 2 e-* 2 / M i 



The single-pole parameters which give stable solutions for A s are c+ = ( 
1.0-22.3 M B ) x 10~ 7 , and for A p are = (-39.8 -19.7 M B ) x 10~ 7 . The results 
are similar to those shown in Fig. 6. 

3 Discussion 

In a recent work by Borasoy and Holstein0 it was found that by including 
resonances in the chiral perturbation theory approach to nonleptonic hyperon 
decays one can obtain fits to data, which does not seem possible if they are 
not explicitly included. In the QCD sum rule method we find it is also 
necessary to include resonances or single pole terms and have done so. 

In summary, we have used QCD sum rules and two 2-point formalisms 
to examine the St and Si nonleptonic decays. We have assumed that pion 
kinetic energy effects are small and have concentrated on the sum rules for 
which we expect these effects to be of order AM/M -C l.We have omitted 
gluonic corrections, in part because we show that the single pole (resonance) 
contributions can be adjusted to fit the data, and the changes in these pa- 
rameters when gluonic corrections are included would not be large. Stable 
and consistent fits to the decay amplitudes are obtained with resonance con- 
tributions of similar magnitude as those used in chiral perturbation theory 
fits. These contributions are not known and are large, so that at the present 
stage the method has no predictive power, which is also true of the chiral 
perturbation theory calculations. Our results, however, parameterize quite 
specific terms in the dispersion relations for the correlation functions that we 
have used, and might be useful for predicting related reactions. 



The work of W-Y.P.H. was supported in part by the National Science 
Council of R.O.C. (NSC88-2112-M002-001Y). The work of E. M. H. was 



Acknowledgements 



15 



supported in part by the U.S. Department of Energy while that of L.S.K. was 
supported in part by the National Science Foundation grant PHY-00070888 
and in part by the U.S. D.O.E. while at the Los Alamos Natonal. Laboratory. 
This work was also supported by the N.S.C. of R.O.C. and the N.S.F. of the 
U.S.A. as a cooperative research project. EMH wishes to thank the Los 
Alamos National Laboratory, the Institute fuer Kernphysik at Juelich and 
the Theory Institute at Tuebingen, where some of this work was carried out 
and completed. In addition, he thanks the Humboldt Foundation for a grant. 
LSK wishes to thank the TQHN group at the University of Maryland and 
the DOE's Institute for Nuclear Theory at the University of Washington for 
hospitality during the period when this work was completed. 



References 

[1] J.F. Donoghue, E. Golowich, and B. R. Holstein, Phys. Rep. 131, 319 
(1986); See also, e.g., J. F. Donoghue, E. Golowich, and B. R. Hol- 
stein, Dynamics of the Standard Model, Cambridge Univ. Press, New 
York (1992). 

[2] P. Zenczykowski, Phys. Rev. D 50, 3285 (1994); D. Horvat, Z. Narancic, 
and D. Tadic, Phys. Rev. D51, 6277 (1995); K. Takamaya and M. Oka 
hep/ph 9811435 (1998) 

[3] B. Barasoy and B.R. Holstein, Phys. Rev. D 59, 094025 (1999) 

[4] M.A. Shifman, A.J. Vainshtein, and V.I. Zhakarov, Nucl. Phys. B147, 
385 (1979). 

[5] B.L. Ioffe and A.V. Smilga, Nucl. Phys. B232, 109 (1984); V.M. Belyaev 
and Ya.I. Kogan, Pisma Zh. Eksp. Teor. Fiz. 37, 611 (1983). 

[6] D.B. Leinweber, Ann. Phys. (NY) 254, 328 (1997). 

[7] E.M. Henley, W-Y.P. Hwang, and L.S. Kisslinger, Phys. Lett. B367, 21 
(1996); addendum Phys. Lett. B 440,449 (1998). 

[8] C.B. Chiu, J. Pasupathy, and S.J. Wilson, Phys. Rev. D32, 1786 (1985). 



16 



[9] L.J. Reinders, H.R. Rubinstein and S. Yazaki, Nucl. Phys. B213, 109 
(1983). 

[10] H. Shiomi and T. Hatsuda, Nucl. Phys. A594, 294 (1995). 
[11] M.C. Birse and B.V. Krippa, Phys. Let. B381, 397 (1996). 
[12] V. Braun and I. Halpern, Phys. Lett. B328, 457 (1994). 
[13] V.M. Belyaev and M.B. Johnson, Phys. Lett. B423, 379 (1998). 



17 



